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ORDINARY-MODE ELECTROMAGNETIC INSTABILITY 
IN COUNTERSTREAMING PLASMAS 
WITH ANISOTROPIC TEMPERATURES 

M. Bomatici* 

Kai Fong Lee 

ABSTRACT 

The instability of the electromagnetic linearly polarized mode 
propagating perpendicular to a uniform magnetic field is studied by 
using the Vlasov equation for a counterstreaming electron plasma 
with anisotropic temperatures. An instability occurs if the stream- 
ing velocity exceeds a certain threshold value which can be below 
that required to excite the electrostatic two-stream instability. It 
is found that temperature perpendicular to the field has a stabiliz- 
ing effect, while parallel temperature enhances the electromagnetic 
instability. Typical growth rates are of the order of the electron 
cyclotron frequency. 


*National Research Council Postdoctoral Research Associate. 
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ORDINARY-MODE ELECTROMAGNETIC INSTABILITY 
IN COUNTERSTREAMING PLASMAS 
WITH ANISOTROPIC TEMPERATURES 

I. INTRODUCTION 

A number of papers have recently appeared in the literature on the linear 

stability of waves propagating perpendicular to the direction of stream motion in 

1-5 

plasmas with relative streaming. In particular, for a system of two identical 
electron plasmas, counterstreaming along an external uniform magnetic field, it 
has been found that in addition to the well-known electrostatic (longitudinal) two- 
stream instability with propagation vector in the direction of streaming, the elec- 
tromagnetic (transverse) linearly polarized mode propagating perpendicular to 

3 4 

the magnetic field can become unstable. ’ The cold plasma theory predicts 
that the threshold streaming velocity required for the excitation of this "modified 
ordinary mode” is cn/a, p , where c is the velocity of light, n and oj p are the elec- 
tron cyclotron and plasma frequencies respectively . The effect of temperature 
perpendicular to the field has also been investigated by means of the fluid equa- 
tions, with the result that the electromagnetic instability is stabilized by per- 

4 

pendicular temperature. 

The purpose of this paper is to study the modified ordinary mode instability 
by taking into account temperatures both perpendicular and parallel to the mag- 
netic field (T I and T n respectively). Our analysis is based on the Vlasov equa- 
tion for counterstreaming plasmas with anisotropic bi-Maxwellian velocity 
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distributions. It is found that an electromagnetic instability can be excited at 
streaming velocities which can be below the threshold value required for the 
electrostatic two-stream instability. It is shown that while T I is stabilizing, in 
qualitative agreement with the result of the macroscopic theory, T| ( enhances the 
instability. This is in marked contrast to the electrostatic two-stream instability, 

g 

which is stabilized by T H but independent of T 1 . 

Typical growth rates for the electromagnetic instability are of the order of 
the electron cyclotron frequency. 

The dispersion relation for the linearly polarized mode for a counterstream- 
ing bi-Maxwellian plasma is given in Section n. Section m presents the instability 
analysis and some graphical illustrations of the results. Finally, in Section IV the 
results are summarized. 

H. DISPERSION RELATION 

Let us consider a homogeneous, infinite plasma in a static, uniform magnetic 
field B 0 . When the plasma is sufficiently hot, the interparticle collisions may be 
neglected and our analysis is based on the Vlasov equation for each charged- 
particle species combined with Maxwell's equations. The system is linearized 
and Fourier transformed in space and time by assuming perturbations propagating 
perpendicular to the magnetic field. We shall assume that the unperturbed veloc- 
ity distribution functions F oj depend only upon the parallel and perpendicular 
components of the particle velocity (V| ( and v A , where subscripts 11 and i refer 
to directions with respect to B 0 ) and are even functions of v„ . In such a case, 
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the general dispersion relation for waves propagating perpendicular to B 0 factors 
into two equations: the dispersion relation for the ordinary mode, which is a 
purely transverse linearly polarized mode, whose electric field is along B 0 , and 
the dispersion relation for the extraordinary mode, which is a mixed transverse- 
longitudinal mode, elliptically polarized in a plane perpendicular to B 0 and con- 
taining the wave vector. We shall restrict ourselves to the investigation of the 

stability of the linearly polarized mode. The corresponding dispersion equation, 

7 

relating the frequency go and the wave-vector *, can be written as 



where 

co . - (477 N . e Vm . ) 1/2 

p j ' ) j/ J / 

and 

= l e jl B o/ m j c : 

Nj is the equilibrium density for particles of type j , J n is a Bessel function of 
order n. 
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We study the above dispersion relation for the case in which particles of the 
same species are counterstreaming along the magnetic field with random (Max- 
wellian) energies in directions both perpendicular and parallel to B 0 . The equi- 
librium velocity distribution F oj is, then, of the form 





( 2 ) 


where u . is the mean directional velocity of particles of species j and 



are their thermal velocities perpendicular and parallel to B 0 respectively (T i i 
and T (| j are non-isotrupie temperatures). The normalization 



holds and no current along the lines of force is associated with a distribution of 
the form (2). From the specific equilibrium velocity distribution (2) and by using 

g 

Weber's second exponential integral, the dispersion relation (1) becomes 




where I is a modified Bessel function and 

n 
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By using the relation 


K )’-" 2 


Equation (3) can be written as 



In the case of no relative motions, u j = 0, Equation (4) yields the dispersion 

9 10 

relation which has been discussed by Hamasaki, and by Davidson and Wu. 

These authors conclude that a purely growing mode exists in an electron plasma 
with /3 ||e > 2 and T le /T |(e < 1, where /3 ))e is the ratio between the electron kinetic 
energy density in the direction of B 0 and the magnetic energy density. In the 
next section, it will appear that in the presence of relative streaming, an instability 
occurs for /3 ||e and T le /T |le which do not necessarily satisfy the above conditions. 

A broadening of the unstable region is as expected since in our system the rela- 
tive streaming motion is a cause of instability as well as the temperature 
anisotropy. 
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in. INSTABILITY ANALYSIS FOR ELECTRON PLASMAS 
A. Instability Criterion 

Henceforth we consider the ion mass to be infinite so that only electron 
dynamics needs to be considered. In such a case an equilibrium distribution 
function of the form (2) describes a symmetrical double beam system, i.e., two 
opposing drifting electron beams which interpenetrate in a static, uniform mag- 
netic field. 11 In order to establish the instability criterion it is convenient to 
use Equation (3), rather than Equation (4), written in the form 



with 


(5) 



where the subscript j = e has been dropped and all quantities refer to electrons. 

g 

Following Hamasaki, the dispersion relation (5) can be analyzed graphically 
by plotting L k (a: 2 ) and R k (<v 2 ) as a function of J 1 for given values of k and look- 
ing for the intersections of the two curves. If R k (0) < L k (0), the smallest root 
for w k 2 ij negative, corresponding to an absolute instability. On making use of 
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Equations (6) and the relation 


CD 

£\<m) = 2 [e“-I 0 GO] . 


the instability criterion can then be expressed as 


(tI +2 ^ _1 )' c2k2 + -p 2 (ti +2 ^) e '“ 


1 "0 T: ■ F ( u 4r **) ; 0 • 


W = 

c 2 n 2 B 2 /87 t 


U = 1 + 2 


^r(us,./i) - * e ~ M I„ (m) • 
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Since F > 0 and 



< 1 , 


two necessary conditions for instability follow from inequality (7): 


<1 + 2 


u 


(9) 


and 


A 


1 +2 


( 10 ) 


Therefore, due to the presence of relative streaming, an instability can occur 
also when T A /r |( > 1 and /3 (| < 2. 

B. Stability-Instability Boundaries 

The equation corresponding to inequality (7), 

T 

1 "iv f Hi • m) * o . (U) 

yields the instability-stability boundaries. Solutions of this equation are obtained 
by looking for values of n for which the curves F(U/3 (| , n) and the straight lines 
(l - 1/U T a /T||), corresponding to different values of the parameters U/3 M and 
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(1/U) T i /T || , are tangent. The dependence of F(U/3|| , upon n is such that the 
curve F(U/3 (| , /x) is tangent to the straight line (l - ( 1/U) T i T |( ) for values of 
im = /j. 0 for which the function F(U/3 n , ij.) has a minimum, i.e. 

dF(U/3 |r ,u)| 


Therefore, Equation (11) is equivalent to the following two coupled equations: 


r 


l - 


T 

ir t: - F ( u ^i’ ^o) = 0 




2 


‘ e ° lot^o) + e 1 i Uo) " 0 


V. 


(12a) 


(12b) 


The above system of equations determines the range of values of the param- 
eters (ll, T A /T (| , /3y ) at the instability boundary. It can be solved numerically to 
express the boundaries of instability in terms of curves separating the stable 
from the unstable region in the T i /T || versus /3 (| plane for various values of U, 
i.e. (u/V H ) 2 . The results are shown in Figures 1 and 2. It appears that while for 
stationary plasmas, u = 0, the minimum unstable value of is 2, in the presence 
of relative streaming, an electromagnetic instability can occur in sufficiently 
low-/3 plasmas. As an example, when u = 2V, ( , the minimum value of ^ for in- 
stability is 0.22. 

It is interesting to recall that the threshold streaming velocity required for 
the electrostatic two-stream (TS) instability is 1.3V (I for a Maxwellian equilibrium 
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g 

distribution function. From Figures 1 and 2 it is seen that it is possible to 

excite the electromagnetic instability alone without simultaneously triggering the 

TS instability. For example, for T a /T h = 1, i.e., an isotropic plasma, and /3 |( ± 3, 

the critical streaming velocity for electromagnetic instability is V, |t which is 

below that required for the TS instability. If T 1 /t || < 1, the threshold is further 

reduced. It therefore follows that conditions can be arranged in the laboratory 

such that the TS instability is suppressed while the electromagnetic instability 

is excited. This conclusion is somewhat different from an earlier one based on 

the fluid equations, in which only the effect of perpendicular temperature was 
4 

considered. 

Corresponding to a set of parameters in the unstable region, there is a 
range of wavenumbers for which the waves are unstable. Such a range is de- 
termined from the instability criterion (7). Since inequality (7) requires the func- 
tion f(U/3h , fi) to be less than a quantity smaller than one, it follows that the un- 
stable range of values of /x is such that /x min < ix </x max . Therefore, an instability 
occurs for wavenumbers k which are greater than a minimum value, k min , and 
smaller than a maximum value, k max . In general, jx min and /z max can be determined 
only numerically. However, by assuming that /x min and ^ max are sufficiently 
smaller and greater than one respectively, the modified Bessel function I 0 (/x), 
which appears in (7) through F(U/3 h , /x), can be approximated with its power series 
and asymptotic expansion, respectively, and expressions for /x min and/x. max can be 
obtained. As an example, let us determine /x min . For /z less than one (/x < 0.2), 
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we can approximate F(U/3|j , /z), at first-order in /z, as 


2/i 


F ( u 4r >*) ' JOT + 1 ' * ■ 


and the instability criterion (7) yields /z > u . , where 


* ii 


1 +2 


IT 


2 


(13) 


Expressions (13) for /x min is consistent with the assumption /z < 0.2 for values of 
parameters jj^Tn , (u/V M ) 2 , /3| ( J such that 


5 < 1 + 2 

* II 


fey- 


2 

7% ’ 


(14) 


Since 


■ i ffi 1 • 


by using Equation (13), we obtain 


k 2 

min 


1 + 2 


u 2 2 


/n \ a 


( 15 ) 


which yields the minimum unstable wavenumber under condition (14). 
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In order to be ablt to relate the results of the present kinetic theory to 
those of the cold theory of the modified ordinary mode, it is useful to solve Equa- 
tions (12) in the limit of low perpendicular temperature, i.e., for m 0 small. In 
this case the exponential and Bessel functions in (12a) and (12b) can be expanded 
in m 0 * At second -order in u 0 we have 


r 


1 l i 

ut: 


* ^0 


= o 


< 


2 

KV 


- l 


3 

+ 1 




2 


= 0 


(16a) 


(16b) 


A quadratic equation in U is obtained by solving one of the above equations with 
respect to m 0 and substituting into the other. Only one of the two solutions of the 
quadratic equation is consistent with the assumption that m 0 be sufficiently small 
than one. The acceptable solution can be expressed as 


u 


2 

c 



+ a‘ 


V, 2 V 2 

v i v || 

“4 T ’ 


(17) 


where 


• [»&)■-]"-{ 


(18) 


a 2 is a positive quantity for 


Vi < 1-5# c 


(19) 
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Equation (17) yields the threshold value for the streaming velocity when the 
temperature perpendicular to the magnetic field is sufficiently low and such that 
condition (19) is satisfied. From Equation (17) it appears explicitly that both the 
magnetic field and the perpendicular temperature have a stabilizing effect, while 
the parallel temperature enhances the instability.* In the cold limit the critical 
value for the streaming velocity becomes ( Q/o. > } c , which is the result derived 
in Reference 3. 

C. Growth Rate 

If the instability criterion (7) is satisfied, the dispersion equation for the 
linearly polarized mode has a negative root for cc 2 , corresponding to a noncon- 
vective instability. By indicating this solution with - "k = Vv > 0 and by taking 
into account only electrons, Equation (4) yields 



In general, the growth rate can be determined by means of a graphical method. 
Simple analytical expressions for the growth rate | y k | can be given only to the 
extent to which one is able to approximate Equation (20). In the approximation 

*Ia Reference 4, by using the fluid equations it was concluded that the temperature stabilizes the 
modified ordinary-mode instability. This result applies only to the particular case in which 
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in which 


1 + 


T 

U T 


\L 

n* 0 


CO 


(nQ) 2 + 7 k 2 




U 


Ti n 2 


+ 7„ 


2e" w I , (fi) 


(21) 


Equation (20) yields 



r 


1_ ut: _f (u/3||, m) 


2e-M, ( M ) 




(22) 


For the above to be a valid expression for the growth rate, the right-hand 
side of (22) has to be real. Hence, in addition to condition (7), which is the in- 
stability criterion for the non-approxi mated Equation (20), the following condition 
has also to be satisfied: 





(23) 


Conditions (7) and (23) are independent conditions. For example, with T L fl |( - 0.8; 
/3 h = 1; (u/Vy ) 2 = 9 and 2 </x <3.2, condition (7), but not (23), is satisfied. In 
Figures 3 and 4 we show some numerical results obtained by using Equation (22). 
The quantity | |/o is plotted versus /x. Ir. agreement with the conclusions of 

Section III B, an instability occurs for /z min < /x < The curves in Figure 3 
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are obtained for (u/V H ) 2 = 2, /3 (| = 2 and four different values of the parameter 
T A /T |t . It appears that the instability is enhanced by increasing parallel tempera- 
ture and stabilized by increasing perpendicular temperature. In Figure 4 curves 
are shown for = 1, T 1 /T 1 | = 0.8 and two different values of the parameter 
(u/V H ) 2 . By comparing these two curves at the same parallel temperature, it 
is seen that the instability is enhanced by increasing streaming velocity. The 
growth rates shown in both Figures 3 and 4 are of the order of the electron 
cyclotron frequency. Since in these examples /3 (| = a> p 2 Vn^Q 2 c 2 ^ 1 and V |( < c, 
it follows that 0 < c j Therefore, the above growth rates for the electromagnetic 

t 

instability are smaller than those corresponding to the electrostatic two-stream 
instability, which are of the order of the electron plasma frequency. 

Instead of (21), the following approximation 


1 + 


n*0 


CO ' 


(nH) 2 +y k 2 


e' 4 I n (M) - 1 


(24) 


has been used by Hamasaki in calculating the growth rate from Equation (20), in 

g 

the ease of no streaming. When condition (24) is satisfied, Equation (20) yields 



00 


p 




(25) 
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The growth rate given by Equation (25) is maximum for values of m = 8U °h 
that the function F(U/3 (| , /n) is minimum in i.e., 



Equations (26) have been solved by Hamasaki in the limit of no streaming, U = 1, 

g 

and with the assumption that n 0 be sufficiently smaller than one. It seems, how- 
ever, that the expression for the growth rate derived in such a way is not con- 
sistent with approximation (24).* We may note, moreover, that approximation 
(24) is unlikely to be satisfied, in general, since it requires one to be smaller 
than the sum of a series whose first terms are of order one. 

IV. CONCLUSION 

In conclusion, we have investigated an electromagnetic instability in a coun- 
terstreaming plasma in a magnetic field and found that conditions exist under 
which the electromagnetic instability is excited independently of the electrostatic 
two-stream instability. We have shown that the effect of temperature perpendicular 
to the field is toward stabilization while parallel temperature enhances the 


•We note also that Hamasalci’s result is a secv .id-order quantity in while io deriving it an ex- 
pansion at first-order in has been used. 
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instability. Due to the presence of streaming, an instability can occur in suf 
ficiently low-/3 plasmas and may then play a role in laboratory experiments. 
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Figure 1. Stability-instability boundaries in the T^/Tu vs /3j ( plane when 
(u/V„)2 = 0; 0.5; 0.75. Unstable regions lie above the curves. 
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